Black Hole Formation in Core-Collapse Supernovae and Time-of-Flight Measurements 

of the Neutrino Masses 

J. F. Beacom* 

Physics Department 161-33, California Institute of Technology, Pasadena, California 91125 
NASA/Fermilab Astrophysics Center, Fermi National Accelerator Laboratory, Batavia, Illinois 60510-0500 

R. N. Boydt 

Departments of Astronomy and Physics, The Ohio State University, Columbus, Ohio 43210 

A. Mezzacappa^ 

Physics Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831 
(October 19, 2000; minor revisions January 8, 2001) 



In large stars that have exhausted their nuclear fuel, the stellar core collapses to a hot and dense 
proto-neutron star that cools by the radiation of neutrinos and antineutrinos of all flavors. Depending 
on its final mass, this may become either a neutron star or a black hole. Black hole formation may 
be triggered by mass accretion or a change in the high-density equation of state. We consider the 
possibility that black hole formation happens when the flux of neutrinos is still measurably high. If 
this occurs, then the neutrino signal from the supernova will be terminated abruptly (the transition 
takes < 0.5 ms). The properties and duration of the signal before the cutoff are important measures 
of both the physics and astrophysics of the cooling proto-neutron star. For the event rates expected 
in present and proposed detectors, the cutoff will generally appear sharp, thus allowing model- 
independent time-of-flight mass tests for the neutrinos after the cutoff. If black hole formation occurs 
relatively early, within a few (~ 1) seconds after core collapse, then the expected luminosities are of 
order Lbh = 10 52 erg/s per flavor. In this case, the neutrino mass sensitivity can be extraordinary. 
For a supernova at a distance D — 10 kpc, SuperKamiokande can detect a i? e mass down to 1.8 eV 
by comparing the arrival times of the high-energy and low-energy neutrinos in v a +p — > e + + n. This 
test will also measure the cutoff time, and will thus allow a mass test of i/ M and v T relative to v c . 
Assuming that and v T are nearly degenerate, as suggested by the atmospheric neutrino results, 
masses down to about 6 eV can be probed with a proposed lead detector of mass Md = 4 kton 
(OMNIS). Remarkably, the neutrino mass sensitivity scales as (D / ' LbhMd) 1 ^ 2 ■ Therefore, direct 
sensitivity to all three neutrino masses in the interesting few-eV range is realistically possible; there 
are no other known techniques that have this capability. 

14.60.Pq, 97.60.Bw, 04.70.-s, 97.60.Lf 



I. INTRODUCTION 

In the past several years, the growing evidence for neu- 
trino oscillations has caused a great deal of excitement 
over the implied nonzero neutrino masses. Oscillation 
phenomena, however, are sensitive only to differences of 
the squared neutrino masses, and thus provide only a 
lower bound on the heavier mass. Without further input, 
the deduced masses can be increased, and the difference 
of masses decreased, providing exactly the same differ- 
ence of squared masses and hence the same oscillation 
phenomena. 

It is therefore of crucial importance to experimentally 
measure or constrain the absolute scale of the neutrino 
masses. Two indirect techniques have been proposed. 
First, the sum of the neutrino masses can be constrained 
by cosmological arguments. The requirement of not over- 
closing the universe gives an upper bound of about 100 
eV jjj . This bound may be improved by considering the 
effects of neutrino masses on the cosmic microwave back- 
ground and the clustering of galaxies; the claimed (fu- 



ture, in some cases) sensitivity is about 1-10 eV M. 
These arguments require that the other cosmological pa- 
rameters are independently known and may not apply if 
the neutrinos decay.Q Second Q), if all of the neutrino 
masses are connected by small measured mass-squared 
differences, then each mass is constrained by the limit on 
the electron neutrino mass from tritium beta decay, now 
about 3 eV || (the direct laboratory limits on the mu 
and tau neutrino masses are 170 keV |(| and 18 MeV [0, 
respectively). If neutrinoless double beta decay were dis- 
covered (i.e., neutrinos were confirmed to have a Majo- 
rana character), then this could anchor the masses at an 



'Furthermore, it has recently been shown that in scenarios 
with a low (MeV-scale) reheating temperature, the neutrinos 
may decouple without reaching equilibrium, leading to a sub- 
stantially lower density than in the usual scenario; this may 
weaken the cosmological neutrino mass bounds by a factor of 
10 or more 
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even lower value ||; the present limit on the combina- 
tion of masses measured in double beta decay is about 0.2 
eV H . Strictly speaking, to use the arguments of Ref. jjj , 
each oscillation signal must first be decisively confirmed, 
including precise measurement of the mixing parameters 
and identification of the oscillated flavors. Until then, we 
must allow for the possibility that there are more relevant 
flavors than there are measured mass-squared differences. 
For example, if the solar neutrino problem is solved by 
v e —* v s oscillations, and if the LSND signal is ruled out, 
then the atmospheric neutrino problem can be solved by 
— ► v T oscillations with a small mass difference and 
large masses, say 10 or 100 eV, as long as 5m 2 ~ 10~ 3 
eV 2 |o). 

Thus, while the indirect constraints on the neutrino 
masses are valuable, it would be much more satisfying to 
have a direct experimental measurement. Presently, the 
best possibility for direct measurement of the mu and 
tau neutrino masses is by time-of-flight differences using 
neutrinos from a Galactic core-collapse supernova. [] At 
lowest order, a neutrino with mass m (in eV) and energy 
E (in MeV) will experience an energy-dependent delay 
(in s) relative to a massless neutrino in traveling over a 
distance D (in 10 kpc): 



At(E) = 0.515 ) D. 



(1) 



The distance is scaled by the approximate distance to 
the Galactic center, though a supernova may be detected 
from anywhere in the Galaxy and its immediate compan- 
ions (e.g., the Magellanic Clouds). SuperKamiokande 
(SK) and the Sudbury Neutrino Observatory (SNO) 
would have good sensitivity to a Galactic supernova, 
collecting of order 10 4 and 10 3 events, respectively (see 
Refs. ]l2|-|l4| and references therein). Unless the decreas- 
ing neutrino luminosity is interrupted by black hole for- 
mation, it should be possible to measure it to very late 
times (some tens of seconds); either outcome would be an 
important probe of the nuclear equation of state jl5],|l6| . 

The primary interest for mass tests is to measure the 
mu and tau neutrino masses relative to the nearly mass- 
less electron neutrino. A neutrino mass test [^2 13 based 
on the average event arrival times (t) can measure a 
mu or tau neutrino mass as small as 45 eV in SK and 



30 eV in SNO. If the mu and tau neutrinos are max- 
imally mixed with nearly degenerate masses, then the 
sensitivity on either mass eigenstate is better by a fac- 
tor of about a/2, i.e., about 30 eV in SK and 20 eV in 
SNO |0,L3[. This test is independent of supernova neu- 
trino emission models, though it does assume that the 
luminosities of the different flavors have similar shapes 
as a function of time, as expected on general grounds 
and also seen in the supernova models |l7| ). In the ab- 
sence of a model, a Kolmogorov-Smirnov test can be 
used to compare the event rates for different flavors of 
neutrinos; it can be shown that this reduces to the (t) 
test Q. Other tests proposed in the literature are ex- 
plicitly model-dependent, and the models have large un- 
certainties. 

While the (t) test could improve the limit on the tau 
neutrino mass by almost six orders of magnitude, it seems 
very difficult to reach the eV range suggested by the 
cosmological and tritium arguments above. It can be 
shown [|l3| that the mass sensitivity generically scales 

with the detector mass Md as 1/M^/ 4 ; therefore, an- 
other order of magnitude in sensitivity in neutrino mass 
would require detectors 10 4 times larger, which seems im- 
possible. It can also be shown (l3) that the sensitivity is 
independent of the distance to the supernova in the case 
where the deduced neutrino mass is compatible with zero 
and only an upper limit is placed. 

In this paper, a comprehensive study that follows our 
recent Letter |1^], we consider the case that the proto- 
neutron star forms a black hole, instead of gradually 
cooling as a stable neutron star. If that happens early 
enough, then the neutrino signals will be abruptly termi- 
nated as the neutrinospheres are enveloped by the event 
horizon of the black hole. In Section II, we discuss the 
conditions required for this to happen and to be observ- 
able, as well as the expected details of the neutrino signal. 
In Section III, we derive the mass effects on the detected 
neutrino event rate in the general case. In Section IV, 
we show how to measure the black hole cutoff time in 
SuperKamiokande, with or without the complicating ef- 
fects of a possible electron neutrino mass. In Section V, 
we show how to make a time-of-flight mass measurement 
of the mu and tau neutrino masses relative to the cutoff 
time measured in SK. Finally, in Section VI, we discuss 
some remaining issues and conclude. 



^As noted by Shrock if neutrinos are mixed, then beta 
decay spectra consist of incoherent contributions from each 
mass eigenstate, where the endpoints depend on the masses, 
and the weighting on the mixing angles. The presence of 
kinks in the spectrum would thus allow direct measurement 
of 77i2 and 77i3 and their mixing angles. In order to experi- 
mentally separate such kinks from an endpoint turnover due 
to mi, the mass differences and the mixing angles must be 
large enough. For light neutrinos, the v e disappearance ex- 
periments presently provide more restrictive limits on these 
parameters. 



II. BLACK HOLE FORMATION AND THE 
SUPERNOVA NEUTRINO SIGNAL 

Before discussing how to measure the neutrino masses, 
we first examine how likely it is that black hole formation 
will truncate the neutrino flux from a Galactic supernova. 
Three questions naturally arise: 

(1) Is the Galactic supernova rate reasonably high? 

(2) Are black holes formed reasonably often in core- 
collapse supernovae? 
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(3) Can black hole formation occur when the neutrino 
fluxes are still high? 

An examination of the evidence reveals that, while the 
uncertainties are large, there is a good chance of sat- 
isfying all three requirements. If so, this could have a 
profound impact on our ability to directly measure all 
three neutrino masses. Before showing how that could 
be done, we address these requirements. 



A. Galactic Supernova Rate 

From studies of other galaxies, we know that about 80- 
90% of supernovae are of the core-collapse type (types 
II, lb, Ic) , which produce a substantial flux of neutri- 
nos po| , pi| . In the following we treat the overall super- 
nova rate without regard to correction for the smaller 
rate of la supernovae. 

A rough estimate of the Galactic supernova rate can 
be made using the historical records. Over the past 1000 
years, 7 Galactic supernovae are known either from his- 
torical records or their remnants pl| , p2| . Probably some 
others in the southern sky were missed because they 
were not visible to or not recorded by the astronomers 
of the time. For example, the recently-discovered su- 
pernova remnant reported in Ref. [^2| is apparently ex- 
tremely close (0.2 kpc) and only about 700 years old, 
but is not found in the historical record. It is therefore 
not unreasonable to estimate that nearby supernovae oc- 
cur at a rate of about 1 /century. Due to obscuration by 
dust, naked-eye supernovae are not visible beyond several 
kpc (the farthest of these 7 was at 4.2 kpc); therefore, 
one must correct for the small fraction of the Galaxy 
surveyed. The Bahcall-Soneira Galactic model [ p3|j24| 
includes somewhat less than 10% of the stars within 
about 4 kpc of Earth; therefore, we estimate the total 
Galactic supernova rate to be about 10/century (see also 
Refs. [|!|26)). 

This estimate of 10/century agrees with the rate given 
by Bahcall and Piran |24|, who make a direct integra- 
tion over the stellar initial mass function, corresponding 
stellar lifetimes, and spatial distribution of stars; their 
calculation is not normalized to the historical rate. It 
also agrees with the nucleosynthesis arguments of Arnett, 
Schramm, and Truran |27|. 

On the other hand, more conservative estimates sug- 
gest that the rate is lower: (3 ± l)/century (|^J2^]. It is 
not clear how to reconcile this with the above estimates 
of 10/century. The estimate based on the historical rate 
and the independent Bahcall-Piran calculation agree, and 
the only element they have in common is the fraction of 
stars nearby. Thus, the most likely fault with these cal- 
culations, if any, is that they assume that the stars that 
explode as supernovae are distributed in the same way 
as other stars. In fact, Refs. |^8|,^9| argue against this 
assumption, and claim that the nearby supernova rate 
is anomalously high due to our occupying a privileged 



position in the Galaxy. 

With coverage over most of the Galaxy over most of 
the past 20 years, no neutrino detectors have reported a 
Galactic supernova ^] (note that SN1987A is excluded 
because it occurred in the Large Magellanic Cloud). 
Taken at face value, this would exclude a Galactic super- 
nova rate of 10/century at about the 85% CL. However, 
an analysis combining all of the experiments has not been 
done, and is needed. A number of these experiments did 
not have full coverage of the Galaxy and/or had signif- 
icant (~ 50%) downtime, and taking this into account 
will yield a weaker constraint. 

LIGO |H| and other novel techniques p2]-p^| may also 
be able to shed some new light on the supernova rate. 

The combined evidence thus suggests a Galactic super- 
nova rate of at least 3/century. 



B. Relative Frequency of Black Hole Formation 

From a theoretical point of view, the relative frequency 
of black hole and neutron star formation (the BH/NS 
ratio) depends on the equation of state of nuclear mat- 
ter |35]] and the supernova mechanism [fj"7f ; further work 
on each is greatly needed. Ideally, appropriate direct ob- 
servational constraints on neutron-star properties could 
be decisive for discriminating between different equations 

of state mm. 

As is well-known, SN1987A in the Large Magellanic 
Cloud (D ~ 50 kpc) was clearly observed by the 
Kamiokande II and 1MB detectors, with 12 and 8 events, 
respectively US The observed duration of SN1987A 
was about 10 s, consistent with a supernova that formed 
a neutron star. No neutron star has been seen yet in 
the remnant, but this may not mean that one is not 
present |j9). Thus if a black hole formed, it evidently 
happened after the neutrino flux died out p0| . The pro- 
genitor mass plays an important role in deciding the ulti- 
mate mass and hence fate of the core. Thus, even though 
SN1987A (progenitor mass ~ 18M ) did not form a black 
hole in the first 10 s after collapse, other supernovae will 
be different. 

Core-collapse supernovae occur only for stars massive 
enough to burn their cores up to iron; this minimum 
mass is estimated to be about 8Mq. It is also gener- 
ally believed that stars above some mass, perhaps 20M Q , 
will always produce black holes instead of neutron stars. 
Bahcall and Piran estimate that supernovae from 
progenitors above 20M Q number about 1/2 of those be- 
low 20M Q . Ratnatunga and van den Bergh ||4l|| , with a 
supernova rate several times smaller, estimate about 1/4 
for this ratio. Fryer |^|| estimates a BH/NS ratio some- 
where between a few percent and 1/4, depending on the 
cutoff progenitor mass; both are strongly affected by the 
uncertainties in the inputs to his supernova code. 

For an assumed stellar initial mass function, predic- 
tions of the remnant mass distribution have been made by 
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Timmes, Woosley, and Weaver p3[ , who find a bimodal 
distribution with peaks at about 1.3M Q and 1.8M©. (In 
some other models, this bimodal distribution is not seen; 
see, e.g., Fig. 3 of Ref. 0]). The bimodal nature in 
this model is due to progenitor masses below and above 
19M Q , which either burn carbon convectively or radia- 
tively, respectively J43| . If the maximum neutron star 
mass is the conventional 2.2M Q Q, then the BH/NS 
ratio ~ Q. However, Brown and Bethe (4^^] argue 
that the maximum neutron star mass is about 1.5Mq, on 
the basis of both an assumed softer equation of state and 
a number of observational constraints. In the Brown and 
Bethe model, progenitors above about 18M Q will form 
black holes, and they independently deduce a BH/NS 
ratio ~ 1. For a maximum neutron star mass of 1.5M©, 
the Timmes, Woosley, and Weaver remnant distribution 
indicates a BH/NS ratio ~ 3 (the upper peak is larger 
than the lower peak) |Q . 

Recent results by Ergma and van den Heuvel E ] 
indicate that the vast majority of progenitors above 
20 — 25M Q produce black holes (this therefore supports 
a much higher BH/NS ratio than the earlier paper of van 
den Heuvel and Habets [Q that suggested a BH/NS ra- 
tio ~ 1/100). This is corroborated by Ref. Q, which 
suggests that the progenitor mass cutoff may be even 
lower. These results thus suggest a high BH/NS ratio. 

It may eventually be possible to address the absolute 
BH/NS ratio observationally via the BH/NS ratio (per- 
haps as large as 10 in a preliminary study [pl[ ) deduced 
from low-mass x-ray binaries, though this also depends 
on the details of the binary evolution. 

Qian, Vogel, and Wasserburg [|2j assumed that the r- 
process production of heavy nuclei occurs in core-collapse 
supernovae and considered the effects of black hole ver- 
sus neutron star formation on the yields. They found 
that the observed r-process distribution may be best ex- 
plained with a very high BH/NS ratio ~ 10. Their results 
require that black hole formation happens early, when the 
neutrino fluxes are relatively high, which will terminate 
part of the r-process production. While their BH/NS ra- 
tio is very large, their hypothesis is supported by recent 
measurements ]53[ ]. Further measurements of r-process 
yields in ultra-mctal-poor stars would be very valuable. 

The accumulated evidence thus supports a relatively 
high BH/NS ratio, so that the next Galactic supernova 
would be likely to form a black hole. 



C. Scenarios for Black Hole Formation 

One scenario for black hole formation in core-collapse 
supernovae occurs if the proto- neutron star mass exceeds 
the maximum neutron star mass. For ordinary neutron- 
rich nuclear matter, this maximum mass is thought to be 
about 2.2M Q pEfl , though there may be significant uncer- 
tainties. This may occur in the initial collapse, or after 
some delay, due to accretion of further mass. The neu- 



trino signal expected in a scenario of this type has been 
studied by Burrows f5~i|| and Mezzacappa and Bruenn |55| l 
(see also a very early paper by Wilson |56|]). In these 
models, neutrino emission was followed until abruptly 
terminated by black hole formation (the results do not 
continue through the short but nonzero black hole forma- 
tion time). Before the cutoffs at 1 - 2 s, the luminosities 
were fairly constant at more than 10 52 erg/s per flavor. 

A second scenario for black hole formation is based on 
a softening of the equation of state in the proto-neutron 
star as the neutrinos are emitted and a phase transition to 
a more exotic state of matter occurs, containing perhaps 
strange mesons or baryons, charged-pion condensates, 
or free quarks. The maximum neutron star mass for 
such exotic nuclear matter is generally lower |^ , ^7| , ^7| , 
perhaps about 1.5Mq. Thus an initially stable proto- 
neutron star may form a black hole after the phase tran- 
sition. The details of the neutrino signal accompanying 
black hole formation in such scenarios have been studied 
by Baumgarte et al. J58| ; see also earlier work fl59|-|6l"[ . A 
detailed study in full general relativity was made of the 
neutrino emission just before and during the formation 
of the black hole. A singularity-avoiding code (6^] was 
used that tracked the emission in the frame of a distant 
observer (i.e., the result is the redshifted, time-dilated 
luminosity that would be seen in a neutrino detector). 
Before the cutoff at about 10 s, the luminosities were 
fairly constant at about 10 51 erg/s per flavor. 

Finally, we discuss two scenarios in which a neutron 
star can become a black hole long after the neutrino flux 
has died away. As such, these scenarios are not of direct 
interest to us. 

In a successful supernova explosion, the outgoing shock 
will pass through the stellar envelope within a few hours 
after core collapse. If a reverse shock forms, it may 
dump matter onto the neutron star and cause it to exceed 
the maximum mass, hence causing black hole formation. 
These fallback scenarios are discussed in Refs. |39| , |63| . 

Gradual accretion onto old neutron stars until the 
maximum mass is exceeded is also possible, and a 
concrete scenario is discussed by Gourgoulhon and 
Haensel Q . Given specific assumptions about the equa- 
tion of state of nuclear matter, they find that in the 
last stages of accretion that the matter will become less 
neutron-rich, and will emit a burst of v e neutrinos with 
(E) ~ 3 MeV. This lasts ~ 0.5 ms until truncated by 
black hole formation. At 10 kpc, we estimate that this 
would cause ~ 3 events above the SK threshold. In fact, 
since their model does not include neutrino opacities, the 
neutrino energy and the luminosity before the cutoff will 
both be lower. Thus, unless the neutron star is very close, 
this would be undetectable (see also Ref. |6j| for a study 
of the sensitivity of LVD). 

Thus there are some concrete models f34[[35 58 in 
which black hole formation occurs early enough to cut 
off the neutrino fluxes when they are still measurably 
high, though the uncertainties are large and depend on 
the details of the supernova models. 
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D. Details of the Neutrino Signal 

In the general case, the observables for each neutrino 
flavor are the luminosity L{t) and temperature T(t) up 
to and during the time of black hole formation. The 
duration of the cutoff must be very short, of order the 
light crossing time 2R/c ~ 0.1 ms. In the most de- 
tailed numerical treatment available j58j, the duration 
of the cutoff is about 0.5 ms. We assume that this will 
be typical for any mechanism of black hole formation. 
For black hole formation at very early times, the ini- 
tial proto-neutron star would be larger than assumed in 
Ref. Ipsj , and one might argue that this would lengthen 
the duration of the cutoff. However, it should be noted 
that what defines the cutoff is the increasing gravitational 
redshift, and this does not become large until the proto- 
neutron star is already very compact. For emission from 
the proto-neutron star, the neutrino gravitational red- 
shifts are moderate; z ~ GM/Rc 2 ~ 0.1. The redshifts 
only become severe during the short cutoff at tBH, when 
z — > oo (using the full expression for z). In any case, 
further modeling of the neutrino signal up to and dur- 
ing black hole formation is needed, especially for black 
hole formation at earlier times. It will be shown below 
that the statistical error in defining the position of the 
cutoff is larger than 0.5 ms; therefore, all of the neutrino 
flavors can be considered to be cut off sharply and simul- 
taneously at a time tsn- These approximations can be 
made because the expected numbers of events during the 
cutoff are less than 1. 

In Ref. HI , some interesting details of the signal dur- 
ing the ~ 0.5 ms cutoff are pointed out. The very last 
neutrinos to be seen will not come from radial paths, but 
rather from unstable circular orbits. It should be noted 
that the calculation of Ref. j5^] only treats neutrinos on 
radial paths. The final decay of the luminosity due to 
the neutrinos on unstable circular orbits is expected to 
be exponential, with a time constant proportional to the 
black hole mass ]66| . Since this time constant is very 
small, r = ZVZGMbh/c 3 ^ 0.04 ms, the number of such 
events (proportional to the disregarded luminosity mul- 
tiplied by this duration) will be negligible for the cases 
considered in this paper. Normally, electron neutrinos 
are emitted from the largest radius and with the lowest 
temperature. At the end of the neutrino signal during 
black hole formation, the electron neutrinos will be cut 
off last and will briefly have a higher temperature than 
the other flavors (due to less gravitational redshift). Un- 
fortunately, all of these details of the transition are not 
observable with the present and proposed detectors, due 
to the limited statistics. For a very close supernova, the 
situation might be different; this will be discussed below. 

The abrupt and simultaneous termination of all fla- 
vors of neutrinos allows a very simple mass test. Since 
the electron neutrino is nearly massless, the termination 
of the v e event rate in SK will signal the black hole forma- 
tion time tsH (the effects of a possible electron neutrino 



mass will be discussed below). Then, any events observed 
after tsn could only have come from neutral-current de- 
tection of time-of- flight delayed, massive v^, v T , v^, and 
v T . We have assumed that the detector background is 
negligible, in the sense that the expected number of back- 
ground events over a typical delay time is <gC 1. 

Before tBH, one would like to measure L(t) and T(t) 
for all of the neutrino flavors. This is straightforward 
for v e and u e , since the detected outgoing lepton carries 
nearly the full neutrino energy in reactions with nuclear 
targets. Since v^, v T , v^, and v T only have enough en- 
ergy to undergo neutral-current interactions, they are in- 
distinguishable. However, for the same reason, they are 
also expected to be produced with the same luminosity 
and temperature. It is not generally possible to measure 
the temperature for these species directly, and it must be 
inferred by the yields on different targets (cross sections 
with different energy dependence sample the spectrum 
differently; see Fig. 3 of Ref. The measurements of 

L(t) and T(t) for the various flavors before tsn, as well 
as the value of tsn itself, are important probes of the 
supernova mechanism and the equation of state p^ , |l6| |. 
They will also be important for measuring the quantities 
needed for the mass measurement, in order to reduce the 
model dependence. 

In the bulk of this paper, we concentrate in the anal- 
ysis on mu and tau neutrino masses near the limit of 
detectability. The mass effects will then not appreciably 
affect the time dependence of the event rate except at the 
sharp cutoff at tsn- In fact, it will be shown that only the 
luminosity and temperature at tBH itself are relevant. In 
the models ]54|,^5],[58| considered, the neutrino luminosi- 
ties and temperatures before tsH are roughly constant 
over the time scales of relevant mass delays. Thus it 
is adequate (and much more convenient analytically) to 
consider that the luminosities and temperatures of all fla- 
vors are constant for some period before the cutoff, and 
that they have simultaneous step- function cutoffs at tBH ■ 
These assumptions will be relaxed below. 

We assume the following temperatures: T = 3.5 MeV 
for u e , T = 5 MeV for P e , and T — 8 MeV for v^, v T , 
Pfj,, and v T [jl7|. This hierarchy is a consequence of the 
different opacities in the proto-neutron star, and the de- 
creasing temperature with increasing radius. The tem- 
peratures in Ref. ]58| were somewhat higher than these 
conventional values, but the authors explain that this is 
probably due to a numerical approximation in the trans- 
port code. 

To illustrate our results quantitatively, we present re- 
sults for two concrete cases. In the first, called "Early," 
black hole formation is assumed to occur a few (~ 1) sec- 
onds after core collapse, when the neutrino luminosities 
are of order 10 52 erg/s per flavor. This case is nomi- 
nally associated with black hole formation by accretion 
onto the proto-neutron star |54|,|5^] . In the second, called 
"Late," black hole formation is assumed to occur within 
several (~ 10) seconds after core collapse, when the neu- 
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trino luminosities are of order 10 51 erg/s per flavor. [ 
This case is nominally associated with black hole forma- 
tion by a softening of the hig h-density equation of state 
in the proto-neutron star | p8| . Direct extraction of the v e 
luminosity from the SN1987A data roughly supports the 
luminosity-time correspondences given here. It should 
be remembered that these are just examples — it will be 
shown that all of the necessary quantities can be mea- 
sured in a realistic situation. 



III. NEUTRINO MASS EFFECTS 

A. Detected Event Rate 

For a constant, normalized, thermal spectrum f(E), 
but a general luminosity L(t), the event rate for neutrinos 
with nonzero mass is 



dN 
~dt 



N T 1 
4ttL> 2 (E) 



dE f(E)a(E)L(t — At(E)) , (2) 



where Nt is the number of targets in the detector, D 
the supernova distance, and (E) the average energy (for 
a Fermi-Dirac spectrum (E) = 3.15T). Generalization 
to a time-dependent spectrum or a shape more general 
than Fermi-Dirac would be straightforward. The argu- 
ment of L(t) is shifted to account for the possible energy- 
dependent delay of a massive neutrino. 

As discussed above, we assume that the luminosity 
and temperature are constant before tsn (for at least 
much longer than the typical delay time), and then van- 
ish abruptly. That is, L(t) = Lbh 0{tsH — t), where 
Lbh is the luminosity at the cutoff. In Eq. (||), we 
need to evaluate this with the delayed argument, i.e., 
L(t - At(E)) = L BH 0{t BH - t + At(E)). For t < t BH , 
the step function is satisfied for all energies, and the event 
rate is 



dN 
~dt 



(t) = C 



L 



BH 



10 51 erg/s 



dEf(E) 



a(E) 



The integral is the thermally-averaged cross section 



dEf(E)a(E). 



(3) 



(4) 



This constant event rate is the same for massless neutri- 
nos; as long as the delays are much less than the total 
duration of the supernova signal and the luminosity is 
constant, then at a given time the number lost by delays 
to later times is compensated by the number gained by 
delays from earlier times. This is not true at the start of 



"'•Recent work of Pons et al. [|L6| suggests that black hole 
formation would occur after a few tens of seconds; however, 
their final luminosities are comparable to what we assume. 



the neutrino signal, but the rise is much less sharp than 
the black hole cutoff, is model-dependent, and is not con- 
sidered further. For t > tsH, there is an upper limit on 
the neutrino energy, which must be small enough for the 
neutrino to be delayed that long after tsn- Then 



dN 
~dt 



(t) = C 



J BH 



10 51 erg/s 



dEf(E) 



a(E) 



(5) 



The upper limit E max is simply the energy that makes 
the argument of the step function 9{t B H — t + At(Ej) 
vanish; using Eq. ([[]), this is 



0.515D 

t — tBH 



(6) 



where the units are as in Eq. ([j]) . Note that the neutrino 
mass and time dependence appear only through the limit 
of integration. If the neutrino energy can be measured, as 
in some charged-current reactions, then the event rates 
for separate ranges of neutrino energy can easily be ob- 
tained. For an H2O detector, the constant C is 



C H 20 = (1.74/s) 



Mr 



1 kton 



10 kpc 
D 



2 r 



1 MeV 



(E) 



(7) 



The constant for a 208 Pb detector can be obtained by 
scaling by the relative number of targets/kton, i.e., 
18/208; therefore 



a 



208Pb 



(0.151/s) 



M 



D 



1 kton 



10 kpc 
D 



1 MeV 



(E) 



(8) 



B. Number Delayed Past tsn 

The expected number of delayed counts Ndei after tsn 
can be determined analytically by integration of Eq. (j^), 
which will be useful when tsn can be measured indepen- 
dently. This is simply 



N de 



f°° , dN.. 
Jt BH dt 



(9) 



C 



L 



BH 



10 51 erg/s 

a{E) 



dt 

tBH "'O 



dE 



'f(E) 



= C 



L 



BH 



10 51 erg/s 



0(t BH - t + At{E)) 



At(E). 



Note that the upper limit on energy in Eq. (Q) was writ- 
ten using the step function 0(tsH — t + At(E)); this step 
function then disappeared in the integration over t. Now 
define 
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where the fa subscript emphasizes that the weighting is 
over f(E)a(E), and not f(E) alone (as for (E)). Then 



N del = (At(E)) Ja 
Lbh 



x C 



10 51 erg/s 



dEf(E) 



a(E) 



(11) 



Recognizing the event rate before (or at) tsn from Eq. [|, 
this becomes 



dN 

Ndel = -ft$BH) X (At{E)) fa . 



(12) 



By use of Eq. (Jlj), we see that Eq. ( pL 0|) simply defines 
the average value of 1/E 2 . By the mean- value theorem 
for integrals, this can be written as l/E%, where E c is a 
constant to be determined. The weighted delay can then 
be expressed as 



(At(E)) 



So 



0.515 



E, 



D. 



The physical significance of the "central" energy E c is 
that it is (to an excellent approximation) simply the 
Gamow peak of the falling thermal spectrum and the 
rising cross section. It can also be determined by numer- 
ical evaluation of Eq. (jl(]). Thus we arrive at the very 
simple and important result: 



dN , 

N d ei = j(1bit) x 0.515 



m , 

e: 1 



(14) 



where the event rate is in s _1 , and the other units are as 
in Eq. (Q). This formula would obviously be true if only 
a single energy contributed and the sharp cutoff in the 
event rate (see Fig. (||) and Fig. ^) below) were simply 
rigidly translated by the delay. But it is remarkable and 
very convenient that it is still true even when there is a 
spectrum of energies and the event rate develops a de- 
caying tail past the cutoff. As derived, this is an exact 
result. 

In the derivation of these results we assumed that the 
luminosity and temperature (and hence also the event 
rate) were nearly constant before tsH, as suggested by 
the results of Ref. [ p^ , ^5|j5^ ] . For an arbitrary event rate, 
a fit can always be made to the event rate before tsH, 
and dN/dt at t bh extracted and used in the formula for 
Ndei- (Below, we also discuss how T and hence E c can 
be extracted from the data). To integrate Ndel as above, 
it is only necessary that the event rate be approximately 
constant over the scale of the small possible mass delays, 
which is a very mild assumption. 

Once the other quantities can be measured, then the 
neutrino mass m is given by Eq. (|l4|). We show how this 
can be done below. 



A. CC Event Rate and Measurement of tsH 

We first consider how well tsn could be measured if 
we knew that m„ e ~ 0. The dominant event rate in SK is 
from v e +p — > e + +n. The cross section p7| as a function 
of the neutrino energy E, including the recoil, weak mag- 
netism, and radiative corrections, is well-approximated 
at typical supernova neutrino energies (where we can dis- 
regard the electron mass) by 



<t(E) = 0.0952 (E - 1.3) 2 (1 - 7E/M) , 



(15) 



for neutrino energies E > 1.8 MeV. In this formula, 
M is the nucleon mass in MeV, and the cross section 
is in 10 -42 cm 2 . For a temperature T — 5 MeV, the 
thermally-averaged cross section (for the sum of the two 
protons in H 2 0) is 44. x 10~ 42 cm 2 . (Thi s is slightly 



smaller than the result used in Refs. 
improved treatment of the corrections 



681, due to an 
Thus, for a 



(13) supernova at 10 kpc as seen in SK (32 kton), the event 



rate due to v e + p 



+ n can be easily calculated. 



Using Eq. (||), the rate just before the cutoff at tsH is 
~ 1500 s _1 in the Early case and ~ 150 s _1 in the Late 
case. After tsn, the rates are zero. We have disregarded 
the 0.5 ms duration |Q of the cutoff, which should con- 
tain about 0.4 events in the Early case and about 0.04 
events in the Late case. Since these are fewer than 1, the 
cutoff can be considered to be sharp. 

How is the cutoff time measured, and what is its error? 
Suppose we have an event rate R(t) measured before the 
unknown cutoff time tsH ■ The time of the last event U as t 
is a lower bound for t bh ■ If t bh were larger than i; ast by 
St, then the number of events expected after ti as t would 
be 5N ~ R(ti as t)St. If Poisson fluctuations caused that 
number 8N to fluctuate to 0, then ti as t would be smaller 
than tBH by about St. This can only occur for SN < 1, or 
St < 1 / R(tiast) ■ Thus, the error in determining the posi- 
tion of a sharp cutoff is generically of the form 1 / R(ti ast ), 
i.e., depending on the number of events N as 1/N. For 
a rate with a tail instead of a sharp cutoff, the error 
in determining the offset time scales instead as 1/ yN ; 
see Ref. |0| for a discussion of the differences. Further- 
more, since ti as t is always less than tsH, a bias correction 
~ 1/R(tiast) should be added to ti ast to estimate tsH- A 
more sophisticated treatment of this problem using order 
statistics |69j7^] yields the same scaling results. 

Thus, for the Early and Late cases, we find that tsn 
will be measured from the charged-current event rate in 
SK with precision slightly better than ~ 1 ms and ~ 10 
ms, respectively. These uncertainties on tsH will have a 
negligible effect on the mu and tau neutrino mass tests 
in the lead detector discussed below. 
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B. Effects of a Nonzero Electron Neutrino Mass 

From the tritium beta decay experiments § , the max- 
imum allowed value of m v (by CPT, the same as mpj 
is about 3 eV. Using Eq. (||), it is straightforward to cal- 
culate the effects of m Ue — 3 eV on the v e + p — > e+ + n 
event rate after tsn m SK. Suppose that tsn were some- 
how known independently. By calculating the event rate 
and integrating, in the Early case we find 21 events af- 
ter the true tsn, with delays as large as about 40 ms. 
// unrecognized, this would bias the extracted tBH to be 
too large, and would seriously degrade the mu and tau 
neutrino mass test (looking ahead to Fig. |J). In the Late 
case, on the other hand, we would have only 2.1 events 
after the true tsn, with delays as large as about 20 ms, 
with less effect on the mu and tau neutrino mass test (see 
Fig-!)- 

However, this potential problem in defining tsH due 
to the unknown electron neutrino mass can easily be 
avoided. In the reaction t' c +j)->e + + )iina Cerenkov 
detector like SK, it is possible to measure the neutrino 
energy by measuring the positron energy and angle. At 
these energies, 



E v ~ {E e + 1.3) 



Ee 

M 



(1 



(16) 



where E e is the positron total energy in MeV, M is the 
nucleon mass, and cos0 is for the positron along the neu- 
trino direction. This follows from the two-body kinemat- 
ics and the small neutron recoil; the full expression is 
given in Ref. p7| . From Eq. (Q), different neutrino ener- 
gies correspond to different delays. At a given time after 
tBH, only energies low enough to have caused a delay that 
large are allowed. The maximum allowed energy, Eq. (|^), 
falls very quickly after tBH, as l/y/t — tsn- Thus, dif- 
ferent ranges of neutrino energy will be terminated at 
different times after tBH, and these can be separated ex- 
perimentally. 

In general, one would use the event rate as a function 
of time and energy, Eq. (||) and Eq. (Jjj), to make an un- 
binned maximum-likelihood fit to the measured neutrino 
energies and times to simultaneously measure both m Ve 
and tBH- However, even without doing that, we can still 
get a good idea of how well we can measure m Ve and t bh 
by splitting the v e +p — > e + +n data into different ranges 
of neutrino energy, which we define as: 

Low : < E < 11.3 MeV, 
Mid : 11.3 < E < 30 MeV, 
High : 30 < E < oo MeV. 

The Low group must be excluded from consideration be- 
cause these events have positron total energy less than 10 
MeV, and can be confused with the 5-10 MeV gammas 
from the neutral-current reaction ^ + 16 O — ► v + 7 + X , 
where X is either n+ 15 or p + 15 N |6£| . In that energy 
range, one would not be able to distinguish delay effects 
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1. The event rate due to u e +p — > e + +n in SK, in the 
Early case, with an assumed distance of 10 kpc. Note that 
only the rate after about tBH is shown, and that the range of 
t — t B H is very short. We took m„ c = 1.8 eV, which is close to 
the minimum mass that can be discerned from this data. The 
labels "Low" (contains 2.4 events past the true tBH), "Mid" 
(4.8 events), "High" (0.5 events), and "All" (7.7 events) refer 
to ranges of neutrino energy defined in the text. 



due to m„ e or and m„ T ■ Generally speaking, m„ c and 
tsH are correlated when extracted from the data (since 
m„ t > has the effect of apparently increasing tsn)- 
However, the High group has much less delay and will 
thus primarily be sensitive to tsn- Then the Mid group 
will principally be sensitive to m„ t , by counting events 
delayed past the tsn determined by the High group. 

In Fig. [I], we show such a possible analysis for the case 
of m Ve = 1.8 eV, in the Early case. The numbers of 
events after the true tsn are: 2.4 (Low), 4.8 (Mid), and 

0. 5 (High). Since in the High group, the number of events 
in the tail is < 1, the cutoff appears sharp and the time 
of the last event (after the bias correction) specifies tBH 
to within the reciprocal of the event rate at the cutoff, 

1. e., about 2 ms. This uncertainty affects the expected 
number in the Mid group by at most ±2 events. Even in 
this case, one can still reliably see a few delayed counts 
after the measured tsn, enough to establish a nonzero 
mass (the statistics are discussed in detail in Section V) . 
We have ignored the 0.4 events expected during the 0.5 
ms duration of the cutoff. 

Thus, in the Early case it will be realistically possible 
to probe electron neutrino masses as small as about 1.8 
eV in SK. The error on the time tsn extracted from the 
same data is not as large as the possible delays (~ 10 
ms, see Fig. |l|), but instead depends on the statistics of 
the High data. Though from the High data alone the 
error on tsn is about 2 ms, we anticipate that a more 
sophisticated fit to all of the data will reduce the error 
somewhat, to about 1 ms. The smallest detectable m„ s 
could probably also be improved slightly. 

In the Late case, the laboratory bound of 3 eV on the 
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electron neutrino mass will generally be stronger than 
that derived from the charged-current signal, and tBH 
will be measured to about 10 ms. 



V. MEASUREMENT OF THE MU AND TAU 
NEUTRINO MASSES 

A. General Framework 

The basic signature of a mu or tau neutrino mass is the 
observation of neutral-current events after tBH- If many 
counts delayed past tsn were observed, then Eq. (|5j) 
could be used to make an unbinned maximum likelihood 
fit to the mass based on how the rate fell off with time. 
The only measurable quantities for any delayed counts 
are their arrival times and their total number, since it is 
not possible to measure the neutrino energy in neutral- 
current interactions. This is simply because not enough 
kinematic variables are measured (the outgoing neutrino 
and the recoiling nucleus are not detected) . In neutrino- 
electron scattering, measurement of the electron energy 
and angle would allow reconstruction of the neutrino en- 
ergy in principle; in practice, the kinematic range of the 
outgoing electron angle is less than the angular resolu- 
tion of the detectors. Thus it is not possible to select 
ranges of neutrino energy as in the measurement. 
While that could be done crudely by exploiting the dif- 
ferent response functions of different targets (see Fig. 3 of 
Ref. |Q ) , it is not necessary if t bh is measured indepen- 
dently in SK. The various neutral-current yields can also 
be used to estimate the v^, v T , z? M , and v T temperature 
T (or, more generally, the spectral shape). 

The test proposed in this paper is to simply count the 
number of neutral-current events after tsH- There is 
a very simple relation between the number of delayed 
counts and the mass, which we quote again because of 
its importance: 



dN ( m 

N del = -j-(t BH ) x 0.515 — 
at \ E r 



D. 



(17) 



where the event rate is in s _1 , and the other units are as 
in Eq. ([!]) . The first important point is that while there is 
a spectrum of energies, only one integral over that spec- 
trum is important, i.e., the one that determines E c . If 
instead we were making a maximum likelihood fit to a 
large number of delayed counts, the precise way the tail 
was filled out would depend on more details of the shape 
of f(E)a(E). The second important point is that after 
consideration of both the supernova neutrino model and 
the detector properties, the only remaining unknown is 
the neutrino mass. The cutoff time tsu can be measured 
in SK. The number of neutral-current counts Ndel will be 
measured between tsH and some suitable stopping point 
that depends on the size of the possible delay effects and 
the detector background rate. The neutral-current event 



rate at tsn due to mu and tau neutrinos will be mea- 
sured with small error since it can be measured over an 
adequately long interval before tsn- As noted, the cen- 
tral energy E c is well-approximated by the Gamow peak 
of the falling spectrum and the rising cross section. Thus 
E c depends on the temperature T; if not assumed from 
theory, this can be estimated from the data, as noted 
above. We assume that the distance D can be deter- 
mined by consideration of the total yield of events or by 
astronomical techniques (although a supernova at more 
than several kpc will be optically obscured by dust, it 
will still be visible at other wavelengths). 

For given measured quantities, the best-fit mass is 



m = E r 



N de 



0.515D %(t BH ) 



(18) 



In the likely case of no delayed events observed, then the 
best- fit mass is obviously m ~ 0. A limit can be placed on 
the mass by considering the largest mass mu m that could 
have faked the massless case. At a chosen confidence 
level, this depends on the largest number of events that 
could have fluctuated down to events. For example, 
using Poisson statistics, an expectation of 2.3 delayed 
counts fluctuates to less than 10% of the time. Then 
mu m is obtained with Eq. ([lj) with Ndel set equal to 2.3. 
If Ndel > is measured, Table | can be used to deduce 
the allowed range of the expected number of counts and 
hence the neutrino mass. Since the fractional error on 
Ndel due to Poisson statistics is large (~ 1/V23 ~ 65%), 
errors on other inputs are expected to be irrelevant. If 
a large number of counts were measured, the Poisson 
relative error would be smaller, and the uncertainties on 
the inputs would play a more important role. 

Dropping all constants of proportionality, we can also 
write mum as 



(E)D 



(19) 



While no longer written in terms of the directly mea- 
sured quantities, this has the advantage of showing the 
dependence on the theoretical inputs more explicitly. For 
a supernova that does not have the sharp cutoff in the 
rate characteristic of black hole formation, the model- 
independent (t) analysis Jl2|;[l3| yields an mu m that is 
independent of the distance D and that scales with the 
detector mass Mr> as l/M^ 4 Jl^]. The different scaling 
with D, and the much more favorable scaling with Afp, 
are consequences of the sharp cutoff in the neutrino flux 
in the present case. 



B. Supernova Neutrino Detection in Lead 

Recently, there has been discussion of building a large 
supernova detector based on 208 Pb |71~|-|74|] . A lead 
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TABLE I. This table shows how a given measured num- 
ber of events N determines a range for the allowed expected 
number of events /i, using Poisson statistics. For the first line, 
fi = 2.3 is the largest expectation that yields N = at least 
10% of the time. For the second line, /j, = 0.1 is the smallest 
expectation that yields N = 1 (or greater) at least 10% of 
the time, and fi = 3.9 is the largest expectation that yields 
N — 1 (or smaller) at least 10% of the time. Successive lines 
are similar. The best-fit fi is shown in parentheses. Using 
Eq. (|l8|), which relates the number of events and the neu- 
trino mass m, the corresponding allowed range in m can be 
determined. Figs. ^ and g can be used for the same purpose. 



measured number 


allowed range of the expected number 


N = 


0.0 < n (~ 0.0) < 2.3 


N = 1 


0.1 < fj, (~ 1.0) < 3.9 


N = 2 


0.5 < m (- 2.0) < 5.3 


N = 3 


1.1 < (j, (~ 3.0) < 6.7 


N = 4 


1.7 < n (~ 4.0) < 8.0 


N = 5 


2.4 < fi (~ 5.0) < 9.3 



detector would observe supernova neutrinos by detect- 
ing neutrons produced through both neutral-current and 
charged-current interactions of the neutrinos with the 
lead nuclei. The neutrons would be produced primarily 
by the neutral-current interactions of u T , D^, and v T , 
because these have the highest temperature. The neu- 
trons could be detected in (for example) a liquid scintil- 
lator doped with ~ 0.1% gadolinium, which has a very 
large neutron-capture cross section, yielding an 8 MeV 
gamma cascade. The neutron capture time in such a 
doped scintillator is very short, of order 0.030 ms J75[ , 
much smaller than the typical mass delays. 

A novel scheme based on a clear solution of lead per- 
chlorate is also being explored |7(|. Neutrons would be 
detected by the 8.6 MeV gamma cascade from capture on 
35 CI, and electrons would be detected by their Cerenkov 
light. 

The neutral-current cross sections for neutrinos and 
antincutrinos on 208 Pb have been calculated by Hargrove 
et al. [[nj and Fuller et al. J77[ . The calculations in this 
paper are based on the Fuller et al. cross section. | While 
the Hargrove et al. and the Fuller et al. results for 
the spectrum-averaged cross sections agree within 20% 
at T = 8 MeV, the underlying calculations are quite dif- 
ferent. As discussed, the cross section uncertainties have 
only a minor effect on the mass test if Ndei is small. 
Nevertheless, a laboratory measurement of the neutrino 



§A very recent calculation by Kolbe and Langanke J78J sug- 
gests a lower neutral-current cross section for neutrinos on 
208 Pb, although the differences with the standard Fuller, Hax- 
ton, and McLaughlin |77| cross section remain unexplained. 



cross sections on lead (perhaps with the ORLAND de- 
tector 1 79 at the Spallation Neutron Source) would be 
valuable. 

Hargrove et al. consider only the allowed contribution. 
The cross section is assumed to be dominated by a narrow 
Ml resonance at 8 MeV, so that 



a(E) ~(E-8 MeV) 2 



(20) 



for neutrino energies E > 8 MeV. However, Fuller et 
al. find that the cross section is dominated by the first- 
forbidden contribution (they also point out some appar- 
ent errors in the Hargrove et al. calculation of the al- 
lowed contribution). Fuller et al. do not provide the 
cross sections as a function of neutrino energy, but in- 
stead only provide thermally-averaged results for various 
assumed spectra. However, it is straightforward to make 
a reasonable fit to a(E) itself. The neutral-current cross 
section is dominated by excitations to the giant dipole 
resonance at 80 MeV/A 1 / 3 ~ 14 MeV. This is just be- 
low the 2-neutron emission threshold, and they find the 
2-neutron emission probability to be very low 5% of 
all neutrons). The cross section can be fit by the form 



a{E) ~(E-14 MeV) 5 



(21) 



for neutrino energies E > 14 MeV. A fit was made to 
the Fuller et al. results, summing the allowed and for- 
bidden (for T = 8 MeV, the latter is about 80% of the 
total) contributions, and summing the results for v and 
v (for either or v T channel). Using this form, the 
leading constant was found to be 2.7 x 10~ 42 cm 2 . Af- 
ter fitting, the thermally-averaged cross sections in the 



were matched to 
lation probability 



first six columns of Table I of Rcf. 
better than 10%. The 1-neutron spa '. 
is approximately independent of energy over the relevant 
range, and can be taken to be 0.90. It should be em- 
phasized that our fits to the cross section and branching 
ratio will only be valid over the limited range of energy 
that we consider. For a temperature T = 8 MeV, the 
thermally-averaged cross section in Eq. (0) for the sum 
of v and v (again, for either or v T ) on 2 Pb, including 
the 1-neutron spallation probability, is about 760. x 10~ 42 
cm 2 . For a supernova at 10 kpc in which the neutrino 
fluxes are not terminated by black hole formation, the 
number of 1-neutron neutral-current events due to v^, 
v T , 9^, and v T , all at T — 8 MeV, is 455 events in 1 
kton of 208 Pb with perfect neutron detection efficiency, 
in agreement with Ref. j77j (who use T — 7.9 MeV). 

It should be noted that the calculations above were 
specifically for 208 Pb, which is 52% of the abundance of 
natural lead. On the basis of the Thomas-Reiche-Kuhn 
sum rule, Fuller et al. |Q argue that the total neutral- 
current neutrino cross sections at these energies should 
scale as a ~ A, where A is the mass number. Thus, the 
total cross sections for the three isotopes of lead should 
be very similar. The position of the giant dipole reso- 
nance changes only as ~ 1/A 1 / 3 , and the 2-neutron emis- 
sion thresholds are 0.7 MeV higher in 206 Pb and 207 Pb; 
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FIG. 2. The results for the combined 1-n neutral-current 
event rate due to u^, v T , and v T in OMNIS. Note that 
only the rate after about tsn is shown. The Early case is 
assumed, with tsn occurring a few (~ 1) seconds after core 
collapse, and luminosities of 10 52 erg/s per flavor at tsn- 
The assumed distance is 10 kpc. Before tBH, there are other 
reactions that produce neutrons; they are not included here, 
and those events will have to be statistically subtracted from 
the measured neutron rate. Maximal f M <-> v T mixing with 
small <5m 2 is assumed, so m ~ m„ 2 ~ m„ 3 . The m = case 
is drawn with a solid line. The m — 6.1 eV case, with 2.3 
events expected in the tail, is the first case that can be reliably 
distinguishable from m = 0, and is drawn with a long-dashed 
line. The results for other masses are drawn with dotted lines. 



therefore, 1-neutron emission will also dominate in these 
isotopes. 

C. Results for a Lead Detector 

In this section, we calculate results for a 208 Pb detector 
that is specified by the number of events expected for a 
supernova at 10 kpc in which the neutrino fluxes are not 
cut off by black hole formation. We assume that the de- 
tector will have ~ 1000 1-neutron neutral-current events 
due to v^, v T , Pp, and v r in this case. A possible design 
for a 4-kton lead detector with about this many events is 
described by Boyd |73|. This design also includes 10 kton 
of iron, with a smaller number of neutral-current events 
(not included in our calculations) . Further refinements in 
the cross section and detector design {73|,[74| (and hence 
the neutron detection efficiency) may affect the mass of 
lead required to meet the design goal of ~ 1000 neutral- 
current events of this type. Using the Fuller et al. |77]] 
cross section, this goal could be met with a 2.2 kton lead 
detector with perfect neutron detection efficiency. We 
refer to this lead detector, whatever its eventual precise 
specifications, as the OMNIS (Observatory for Multifla- 
vor Neutrinos from Supernovae) detector. 
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FIG. 3. As in Fig. except that the Late case is assumed, 
with tsn occurring within several (~ 10) seconds after core 
collapse, and luminosities of 10 51 erg/s per flavor at tBH- The 
m — 19.2 eV case, with 2.3 events expected in the tail, is the 
first case that can be reliably distinguishable from m — 0, and 
is drawn with a long-dashed line. Note the changes of scale 
on the axes. 



In the following, we assume a supernova distance of 10 
kpc. Using the product of the thermally-averaged cross 
section and the branching ratio given above, the event 



rates due to neutral-current detection of u^, v Tl 



Vp, and 



v r can easily be calculated with Eq. (0) and Eq. (||). 
These rates are shown in Fig. || and Fig^| for the Early 
and Late cases. Recall that the luminosities and cutoff 
times chosen are simply examples; in a real case, the 
relevant quantities will be measured, not assumed. In 
particular, tsn will be measured using the D e +p — > e + +n 
events in SK. 

In Fig. H for the Early case and in Fig. [H] for the Late 
case, the number of delayed events Ndei (that is, i^, v Tl 
D^, and v T events after tsu) is shown versus the neutrino 
mass. The points are from direct numerical integration 
of Eq. ( H) , and the solid line is the simple analytic result 
of Eq. (|l7|). Note that E c = 40.7 MeV is calculated using 
the Gamow peak of f(E )a( E), and is not fitted. 

In order to use Eq. (|18| ), a minor correction to the 
measured event rate before tsn must be made. In a lead 
detector, one expects to measure just the total neutron 
rate. Thus the expected contributions from the charged- 
current 1-neutron and 2-neutron events will have to be 
statistically subtracted, along with the contributions of 
v e and D e to the neutral-current rate. The subtracted rate 
of neutrons before tBH is about 20% of the total \!% . 

The cross section normalization appears only in the 
event rate, where it is multiplied by Lbh, which is a 
priori unknown. Only their product, in the form of the 
measured event rate, is needed in Eq. ([l8|). The cross 
section shape only affects E c . Using a Fermi-Dirac spec- 
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FIG. 4. The expected number of delayed counts Ndei (those 
after tsn , due to the mass effects) in OMNIS as a function of 
the neutrino mass. The calculation uses the same assumptions 
as in Fig. g, the Early case. The points are obtained by 
direct numerical integration. The "+" indicates the smallest 
discernible mass at the 90% CL. The solid line is obtained 
with Eq. @, using E c = 40.7 MeV, the Gamow peak energy. 



trum with temperature T = 8 MeV, then E c ~ 41 MeV 
using the Fuller et al. |77j cross section given above, and 
E c ~ 35 MeV using the Hargrove et al. |7lJ cross section 
given above; this is a negligible difference. The spectral 
temperature T (nominally 8 MeV) of the mu and tau 
neutrinos at the time of the cutoff is a priori unknown, 
perhaps by ±25%, and this also affects E c . The heavy- 
flavor temperature can be estimated from the data by the 
yields on different targets (see Fig. 3 of Ref. |14|| ), and 
this may reduce the uncertainty on T. Thus, in terms 
of impact on the measurement of the neutrino mass, the 
uncertainties in the thermally-averaged neutral-current 
cross section on 208 Pb are of less importance than the 
Poisson counting error. 

Using Figs. ^ and [|, we obtain mass sensitivity as low 
as 6.1 eV in the Early case and 19.2 eV in the Late case. 
These are the first masses that can be reliably discerned 
(90% CL) from the massless case, since they correspond 
to at least 2.3 expected events after tsH- Larger masses 
give even more delayed events, and hence are easier to 
measure. In these results, we have assumed that and 
v T are maximally mixed, with Sm 2 ~ 10 -3 eV 2 , as sug- 
gested by the atmospheric neutrino results [10|, so that 
both contribute to Ndei- The results for the neutrino 
mass will then apply to the two relevant mass eigen- 
states. If we do not consider this mixing, then perhaps 
only the tau (or mu) neutrinos will have a mass and be 
delayed. Then Ndei is half as large as assumed here, and 
by Eq. (|l8|), mu m is \/2 larger. Since assuming that only 
one neutrino is massive is the most conservative possibil- 
ity, the deduced limit would in fact apply for either of 
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FIG. 5. As in Fig 
assumptions of Fig. t 
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^, but for the Late case, and with the 
Note the change in the horizontal scale. 



the mu and tau neutrino masses. 

Finally, we discuss some sources of error for the number 
of delayed events Ndei m a 208 Pb detector, all of which 
are negligible. We ignore possible detector backgrounds 
over the short time scale of possible delays. The duration 
of the cutoff is about 0.5 ms |5S[| ; taking that into account 
would make N de i larger by ~ 0.5 x 200 x 0.0005 = 0.05 
events in the Early case and 0.005 events in the Late case. 
As noted, the uncertainty on tBH from SK is assumed 
to be about 1 ms in the Early case and 10 ms in the 
Late case. From Fig. (||) and Fig. (|J), this uncertainty 
can be seen to change the expected number Ndei by ~ 
±200x0.001 = ±0.2 events in the Early case and ~ ±20x 
0.010 = ±0.2 events in the Late case. Even with m Ve < 
1.8 eV determined in SK, there can still be some v e and 
P e events (charged- and neutral-current on 208 Pb) after 
the true tBH- In the worst case, assuming no tagging 
on 2-neutron events or events with an electron, the v e 
and v e events contribute about 20% of the total neutron 
rate before tsn- Assuming m Ve — 1.8 eV and E c ~ 30 
MeV, then the number of these events after the true tBH 
is ~ 50 x 0.515 x (1.8/30) 2 = 0.09 in the Early case 
and 0.009 in the Late case. For a larger lead detector 
or a closer supernova, some of these errors could become 
relevant. 



D. Results for SNO 

The principal neutral-current reactions available in 
SNO are v + d —* v + p + n and D + d — > v + p + n, de- 
tected by neutron capture. For a supernova at 10 kpc in 
which the neutrino fluxes are not truncated by black hole 
formation, 485 events are expected, of which 400 would 
be caused by u^, u T , v^, and v T |R3]. Perfect neutron 
detection efficiency is assumed. Before tsH, the neutral- 
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current event rate due to these flavors may be obtained 
by scaling the 208 Pb results by 400/1000, the ratio of the 
total numbers of events expected for a supernova that 
does not form a black hole. This works simply because 
both the event rate before tsH and the total number of 
events have the same dependence on a e f f and the number 
of targets. Then, using Eq. (|J) with E c = 32 MeV @ 
and Ndei = 2.3, we obtain mu rn = 8 eV in the Early case 
and vriura = 24 eV in the Late case. 

However, it may not be possible to reach this sensi- 
tivity in practice due to the long neutron capture time 
in heavy water (an exponential distribution with time 
constant t„). The value of t„ depends on the neutron 
capture technique: with the dissolved MgCl2 salt, r„ ~ 4 
ms; with the 3 He counters, r„ ~ 16 ms; and with pure 
D2O, r n ~ 35 ms pel . The effect of this smearing is to 
delay events after tgn even in the massless case: 

dN 

N de i -> N del + -7-(t BH ) x r„ . (22) 
at 

For the Early case, this adds 0.8(r„/10 ms) events after 
tsH- Thus, unless the salt is used, the neutrino mass 
sensitivity of SNO will be degraded because events after 
tBH can be delayed by either and v r mass effects or 
the nonzero neutron capture time. 



E. Results for SK 

The first set of neutral-current reactions available in 
SK are those on 16 discussed above that yield a 5-10 
MeV gamma in the final state ||68[ . For v^, v T , v^, and 
D T , 710 events in total are expected for a supernova at 10 
kpc [l2] ] . Before tBH , the neutral-current event rate may 
be obtained by scaling the 208 Pb results by 710/1000. 
In practice, this event rate will be obtained from the 
measured one by statistically subtracting the comparable 
rate due to low-energy v e + p — > e + + n events, which are 
indistinguishable in SK. Using Eq. ( |l8|) with E c = 60 
MeV |L2| and Ng e i — 2.3, we obtain mu m = 11 eV in the 
Early case and mn m = 34 eV in the Late case. 

However, it may not be possible to reach this sensitiv- 
ity in practice because of the low-energy D e + p — > e + + n 
events after tsH, of which there can be as many as 2.4 
in the Early case, due to the limited sensitivity to m„ e 
in SK. Furthermore, the very steep cross section on 16 
is much more sensitive to the temperature or the spec- 
tral shape in general (see Fig. 3 of Ref . jl4j ) , and so this 
result is more model-dependent. Thus the mu and tau 
neutrino mass sensitivity of SK using the neutral-current 
reactions on 16 will be limited. 

The second set of neutral-current reactions available 
in SK are v + e~ — * v + e~ and v + e~ — * D + e~, for 
which 120 events due to v T , and v r are expected 
for a supernova at 10 kpc [Q. Before tBH, the event 
rate for these reactions may be obtained by scaling the 
208 Pb results by 120/1000. One must first subtract from 



the measured event rate events due to v e + e~ — > v e + e~ , 
v e + e~ — > v e + e~ , and v e + p — ► e + + n in the forward 
cone. The unwanted events dominate the signal before 
tBH by a factor of ~ 5, so the statistical subtraction will 
introduce some error. If this effect can be ignored, then 
using Eq. © with E c = 25 MeV [§J| and N del = 2.3, 
we obtain mu m = 11 cV in the Early case and mu m = 34 
eV in the Late case. 

However, it may not be possible to reach this sensitiv- 
ity in practice, again because of the limited sensitivity 
to m„ c , which can allow otherwise indistinguishable v e 
and v e events after tsn- In the Early case, we estimate 
that there could be ~ 0.9 such events after tsH- Thus, 
the mu and tau neutrino mass sensitivity of SK using the 
neutral-current reactions on electrons will also be limited. 



VI. DISCUSSION AND CONCLUSIONS 

A. Distance Dependence of the Neutrino Mass 
Sensitivity 

Throughout this paper, we have assumed that the next 
Galactic supernova will be at a distance of 10 kpc. In 
the Bahcall-Soneira Galactic model 25%, 50%, 

and 75% of supernovae are within about 7, 10, and 14 
kpc of Earth, respectively. If the events during the short 
(~ 0.5 ms) cutoff can be disregarded, then the results 
for other distances can be scaled with Eq. (p"9j), and are 
shown in Fig. ^| and Fig. |?]. Other errors, for example the 
error in N^ei that comes from the small error on tsH> are 
independent of D in their relative importance. 

A close supernova at 1 kpc would obviously have 100 
times as many events as we have assumed, and would 
naively have mass sensitivity about 3 times better than 
at 10 kpc, i.e., about 2 eV in the Early case. However, 
there could be a number of events during the short cut- 
off that would make defining tsn more difficult than 
for a more distant supernova (even assuming that the 
high event rate in SK does not saturate the detector). 
Assuming a ~ 0.5 ms duration Q, there could be 40 
such events in SK in the Early case and about 4 in the 
Late case. Note that these are estimated simply by the 
area of the triangle with height given by the event rate 
at tsH arid width given by 0.5 ms. In fact, the neu- 
trino temperatures are falling rapidly during these 0.5 
ms, due to increasing gravitational redshift; taking that 
and the detection threshold into account would reduce 
these numbers. Even if tsn could be defined with neg- 
ligible error, there could still be neutral-current events 
after tBH due to the ~ 0.5 ms duration of the cutoff: 
perhaps 5 events in OMNIS in the Early case and 0.5 
events in the Late case. Again, these are conservatively 
large estimates. The presence of events during the cut- 
off would weaken the mass sensitivity, and it would no 
longer decrease with decreasing distance. However, the 
real behavior of the luminosity and temperature during 
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FIG. 6. The mass sensitivity as a function of the super- 
nova distance (solid lines) , for the Early case, for 
measured in the OMNIS detector, and for m„ c measured in 
SK. This figure is appropriate if Ndel — is measured and 
only a limit is being placed on the neutrino mass (if Ndel > 
is measured and hence a nonzero mass is discovered, see Ta- 
ble I). The dashed line is the present laboratory upper limit 
on m Ve fj . In using Eq. ( |l9| ) to make this figure, we assumed 
that the events in the ~ 0.5 ms tail can be disregarded. De- 
pending on the unknown details of the tail, this assumption 
will break down at perhaps ~ 3 kpc and the mass sensitivity 
will not improve further with decreasing distance. 



the cutoff is not well known, and further modeling along 
the lines of Baumgarte et al. is needed. 

For an extremely close (and hence rare) supernova, 
e.g., Betelgeuse at ~ 0.1 kpc, the possibilities are even 
greater, particularly for exploring the process of black 
hole formation |38|,[59| , provided that the neutrino obser- 
vatories can accommodate the enormous event rates. 



B. Neutrino Oscillations 

While a full discussion of neutrino oscillations is be- 
yond the scope of this paper, we make a few brief com- 
ments. Oscillations of v T are not important in the 
sense that these flavors cannot be distinguished experi- 
mentally. The atmospheric neutrino results suggest that 
both are massive, with a small mass difference and a large 
mixing angle |lCf| ; if so, the measured mass corresponds 
to the nearly degenerate mass eigenstates. Oscillations of 
Vhi v t ~ > v s will decrease the number of neutral-current 
events; this is irrelevant in the sense that the mass mea- 
surement depends on the measured, not predicted, event 
rate at tBH- Oscillations of v^v T <-» v e (and their an- 
tiparticles) can in principle complicate the mass tests. 
However, because of the higher temperature for and 
v T , such oscillations would greatly increase the number 
of charged-current events and would harden the electron 
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D [10 kpc] 

FIG. 7. As in Fig. ^, but for the Late case. Because of the 
lower luminosity, the mass sensitivity may flatten out only 
below about ~ 1 kpc. Note the change in the vertical scale. 



or positron spectrum; see, e.g., Refs. [|77|]81| 1. If evidence 
of such oscillations were seen, the formalism presented 
here could easily be enlarged to include oscillations. The 
positron spectrum from v e + p — + 
appears to exclude large Vp,,v T <-* D t 



+ n from SN1987A 
mixing |82|. 



C. Conclusions 

If a black hole forms early in a core-collapse supernova, 
then the fluxes of the various flavors of neutrinos will be 
abruptly and simultaneously terminated when the neu- 
trinospheres are enveloped by the event horizon. For a 
massive neutrino, the cutoff in the arrival time will be 
delayed by At ~ (m/E) 2 relative to a massless neutrino. 
The SK detector can measure both tsn and m„ E by the 
arrival times of low- and high-energy v e + p — > e + + n 
events, for which the neutrino energies can be measured. 
The mu and tau neutrinos are detectable only by their 
neutral-current interactions, in which their energies are 
not measured. However, their masses can be measured 
by counting the number of these neutral-current events 
detected after tBH- 

The mass sensitivity depends on the supernova neu- 
trino luminosity Lbh at cutoff, the distance D, and the 
detector used. For luminosities of 10 52 erg/s per flavor at 
cutoff (the Early case), and a distance of 10 kpc, SK will 
be able to measure an electron neutrino mass as small as 
1.8 eV and OMNIS would be able to measure m v ~ m„ T 
as small as about 6 eV. These results are perhaps even 



slightly conservative, as the luminosities in Refs. |54 55 
were in fact a few times larger than assumed in the Early 
case. As discussed, the mu and tau neutrino masses were 
assumed to be degenerate because of the atmospheric 
neutrino results Jig] ; in this case the masses are really 
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those of the relevant mass eigenstates. 

Using the neutral-current channels in SNO and SK, 
the neutrino mass sensitivity is nominally ~ 10 eV for 
each. However, it appears that various practical effects 
will degrade those results. 

For other luminosities, distances, and detector masses, 
the mass sensitivity scales as in Eq. (|l9|), i.e., 



mii 



D 



LbhMd 



(23) 



This should be contrasted with the case in which the neu- 
trino luminosities are not truncated by black hole forma- 
tion, where mn m ~ 1/M^/ 4 and is independent of D. 

As we have discussed, there seems to be a good chance 
that the ongoing and proposed neutrino detectors can 
observe the truncation of the neutrino signals caused by 
black hole formation in a Galactic core-collapse super- 
nova. This would have profound consequences, even if 
no delayed events were observed and only limits were 
placed on the neutrino masses. Besides the obvious as- 
trophysical importance of such an observation, this could 
improve the limit on the tau neutrino mass by a factor of 
almost 10 7 . Moreover, the technique discussed in this pa- 
per is the only known possibility for direct measurement 
of the Vn and v T masses (either Dirac or Majorana) in 
the crucial eV range suggested by the indirect neutrino 
mass tests [CPJ3l discussed in the Introduction. 
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